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ABSTRACT 

We study the scaling relations between the luminosities, sizes, and rotation velocities of disk galax- 
ies in the SFI++, with a focus on the size-luminosity (RL) and size-rotation velocity (RV) relations. 
Using isophotal radii instead of disk scale-lengths as a size indicator, we find relations that are sig- 
nificantly tighter than previously reported: the correlation coefficients of the template RL and RV 
relations are r = 0.97 and r = 0.85, which rival that of the more widely studied LV (Tully-Fisher) 
relation. The scatter in the SFI++ RL relation is 2.5-4 times smaller than previously reported for 
various samples, which we attribute to the reliability of isophotal radii relative to disk scale-lengths. 
After carefully accounting for all measurement errors, our scaling relation error budgets arc consis- 
tent with a constant intrinsic scatter in the LV and RV relations for velocity widths logVF > 2.4, 
with evidence for increasing intrinsic scatter below this threshold. The scatter in the RL relation is 
consistent with constant intrinsic scatter that is biased by incompleteness at the low-L end. Possible 
applications of the unprecedentedly tight SFI++ RV and RL relations are investigated. Just like 
the Tully-Fisher relation, the RV relation can be used as a distance indicator: we derive distances 
to galaxies with primary Cepheid distances that are accurate to 25%, and reverse the problem to 
measure a Hubble constant Hq = 72 ± 7 kms -1 Mpc" 1 . Combining the small intrinsic scatter of our 
RL relation (e^ = 0.034 ± 0.001 Iog[ft — 1 kpc]) with a simple model for disk galaxy formation, we find 
an upper limit on the range of disk spin parameters that is a factor of ~ 7 smaller than that of the 
halo spin parameters predicted by cosmological simulations. This likely implies that the halos hosting 
Sc galaxies have a much narrower distribution of spin parameters than previously thought. 
Subject headings: galaxies: spirals - galaxies: fundamental parameters - galaxies: photometry - 
cosmological parameters 



1. INTRODUCTION 

The observed correlations between spiral galaxy lumi- 
nosities (L), sizes (R) and rotation velocities (V) have 
long been exploited to provide insight into the large- 
scale galaxy distribution and the natur e of disk galax- 
ies th emselves. Calibrated Tully-Fisher (|Tullv fc Fisherl 
|1977|) scaling relations built from the distance-dependent 
L and the distance-independent V (we henceforth refer 
to this as the LV relation) have been extensively used 
as distance indicators, while scatter in L caused by de- 
viations from the Universal expansion probe the large- 
scale mass distribution. The tightness of the LV rela- 
tion as well as that between R and V (the RV relation) 
and between R and L (the RL relation) over decades 
in mass also provide important clues to the process of 
galaxy formation. Most theories connect the structure 
of galact ic disks to that of their parent dark matter ha- 
los (e.g. iFall fe Efstathioul 119801 see below), themselves 
well-constrained by simulations of halo assembly. A care- 
ful consideration of the scaling relations in the context 
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of galaxy formation thus probes how baryons collapse to 
form disks as well as the properties of the dark matter 
halos that host disk galaxies. 

The LV relation has been used routinely over the 
last 30 years to compute distances to nearby galax- 
ies and their peculiar velocities relative to the Hub- 
ble flow. With typical accuracies of 15 — 20%, LV 
distances have been used successfully for example to 
determine the 3D structure of the Virgo cluster (e.g. 



Tully & Shaval [19841: IPierce fc Tullvl U98a lYasuda etatl 
1997t iGavazzi et al.lll999t ) and to compute the value of 



Hq, with accuracy comparable to cosmology experiment- 
derived values (e.g.lGiovanelli et al.l[l9 97a; Sakai fc et al.1 



120001: iMasters et al.ll2006l . hereafter M06). Because of the 
usefulness, accuracy and relative ease of measurement 
of LV distances, large galaxy samples have been assem- 
bled for this pur pose at most optical and near-infrared 
rave l engths (e.g. | Aaronson et al.ll l982: M athewson et al.1 



[1991 IWiilick et al.H1997t ITheureau et al.lll998D . 

The SFI++ (jSpringob et all 120071 hereafter S07) is a 
sample of ~ 5000 local galaxies, one of the largest of its 
kind, and represents the culmination of 15 years worth of 
work by Giovanelli, Haynes and collaborators (hereafter 
collectively referred to as "the Cornell group") towards 
measuring peculiar velocities using LV distances. The 
SFI++ builds on the earlier Spiral Field I-band (SFI: 
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IGiovanelli et"aI1[T99l fl995t IHavnes et alJH999alfbh . Spi 



ral Cluster I-band (SCI: IGiovanelli et al 1ll997bn . and Spi- 
ral Cluster I-band 2 fSC2: iDale et al.lll999a flbT) samples. 
It also includes additional data from the previously un- 
published Spiral Field I-band 2 sampl e (SF2) and from 
the so uthe rn hemisphere surveys of iMathewson et aTl 
(|1992t) and IMathewson fe Forcfl (jl996l) . with photometry 
re-processed to match that of the other catalogs. The 
sample contains both field and cluster galaxies, mostly 
of type Sc. A subset of the cluster galaxies comprise 
a template sample used to calibrate the LV relation 
via th e "basket-of-clusters" technique (|Giovanelli et al.l 
Il997bl M06). The remaining catalog members comprise 
the nontemplate sample, from which peculiar velocities 
are estimated by application of the template LV relation 
(S07). 

The LV, RV, and RL scaling relations are also impor- 
tant diagnostics of galaxy formation: L is fundamentally 
a tracer of stellar mass, R a tracer of disk specific an- 
gular momentum, and V a tracer of (dark + luminous) 
mass. In the standard model for galaxy formation, 
disks form through the dissipative collapse of hot gas 
within the potential wells of their dark matter halos, 
conse r ving specific angu l ar momentum dWhite fe Reesl 
l97l iFal fc Efstathioul [19801 iDalcanton et al.l 119971; 
Mo et al.l I1998D . This scenario relates L, R and V 



of galactic disks directly to the virial properties of 
their parent halos, which can b e reliably measure d 
from collisionlcss simulations (e.g. iMaccio et al.ll2008h . 
These models have been extensively applied to the LV 
relations origin a lly derived as di s tance indicators (e.g. 
Mo et all 119981: IDalcanton et all 119971: Ivan den Boschl 



19981: ISomervilie fe Primackl 119991: Ivan den Boschl 120001: 



Navarro fe Steinmetzl I200C ; iFirmani fe Avila-Reesd 
and to the RL relation to a lesser extent 
(iSaloeter fc Hoffman! [19961 Ide Jong fc Lacevl l2000l: 



Graham! 120021: iShen et all 120031: fAvila- Reese et ail 
20081) . Thev can typically be "tuned" to match the 
slopes, zero-points and scatters of a given relation. 
However, simultaneously reproducing all scaling rela- 
tions and the galaxy luminosity function is a notoriously 
difficult problem, as it involves both cosmological 
initial con ditions and the physic s of star formation and 
feedback (Ivan d e n Boschl 120001: ISomervilie fc Primack 



li M iBell et al.l 120031: iCroton et al.l 120061 iDutton et al 
2007L e.g.). 

In light of their importance in constraining galaxy 
formation models, scaling relations are now being 
constructed spec ifi cally for this applicatio n (e.g . 
Kannappan et al. 120021: iKauffmann et al.l 120031: 
Shen et all l2003t iZavala et all 120031: IPizagno et al I 



20051: ICourteau" et al .1 12007t IPizagno et all I2007TT Of 

most relevance to the present study, ICourteau et all 
(|2007l hereafter C07) and IPizagno et all (|2005l . hereafter 
P05) have constructed the LRV scaling relations from 
photometry and kinematics comparable to that in 
the SFI++. C07 combines subsets of four extant LV 
distance samples into a catalog of 1300 galaxies with V 
estimated from optical rotation curves, R and L derived 
from /-band scale-lengths and apparent magnitudes 
respectively, as well as as an estimate of R from 2MASS 
.ftf-band effective radii. The raw magnitudes, velocity 
widths and scale lengths are inhomogencous, but C07 
apply the same inclination and extinction corrections 



to all galaxies. By contrast, P05 use a significantly 
smaller sample of 81 disk-dominated galaxies, but with 
homogeneously measured V from op tical rotation curves, 
L from SDSS (jYork fc et alJl2000f ) i-band magnitudes, 
and R from i-band scale-lengths via 2-dimensional 
bulge-disk decompositions. Additionally, they use SDSS 
g — r colors to estimat e stell ar masses. 

IDutton et ail (|2007l) and IGnedin et all (|2007f ) apply 
galaxy formation models to simultaneously fit the C07 
and P05 LRV relations, respectively. Both claim to 
match the scaling relations, their residuals and scatter 
as well as the galaxy luminosity function, but not with- 
out invoking n on-standard parame ters such as low mass- 
to-light ratios (IGnedin et all 120071) or processes such as 
halo expansion (|Dutton et alll2007D . These studies have 
clearly demonstrated the need for large, homogeneous 
samples of galaxies from which scaling relations with 
well-understood residuals, scatters and uncertainties can 
be derived. 

The exponential scale-length of galaxy disks is 
typically adopted as the measure of R in disk galaxy 
scaling relations. However, disk scale lengths are no- 
toriously difficult quantities to measure. For low- 
inclination systems (i < 50°), tests of measured scale- 
length rel iability typically return uncertainties of 10- 
20% fe.g. iSchombert fc Bothu^ll987t |B~vun fc Freeman 



19951: Ide Jond 119961: iMacArthur et alll2003l: iFathi et al 



20101 ) , while comparisons of reported by different au- 
thors for the same galaxy reveal a scatter of ~ 25% 
(|Knapen fc van der KruTtlll991l : lMollenhoflll2004f ). Most 
scaling relation samples contain galaxies with i > 60° 
to mitigate uncertainties in the inclination correction 
required to estimate V. In the optical, disks are 
parti ally opaque in the region where is measured 
(e.g. IGiovanelli et all ll994. 1995), and internal extinc- 
tion is a significant additional source of uncertainty in 
the derivation of R in the scaling relations. R ecipes 
to co r rect rg for intern a l extinction exis t (e. 
I99l IBvun et all fl99i IGiovanelli et"ail [199 



Bvunl 



1995 



Graham 2001 ; Masters et al. 120031: IMollenhoff et all 
20061: iGraham fc Worlevl l2008f ). and imply that large 



(> 20%), uncertain (by ~ 30%) corrections are required 
even at moderate inclinations of i ~ 70°. 

Given the large and potentially systematic uncertain- 
ties in r^, it is perhaps not surprising that the scaling 
relations constructed using this parameter show weaker 
correlations and larger scatters than the LV relation 
(e.g. C07, P05), and have been correspondingly less well- 
studied than the latter in both the distance indicator and 
galaxy formation contexts. It therefore seems worthwhile 
to explore other measures of R than r^, such as isopho- 
tal radii, to const ruct the LRV s calin g relations: pre- 
limina ry work by iSpekkensl (|2005l ) and iSaintonge et all 
(20081) has demonstrated promise in this approach. 

While originally designed to compute peculiar veloci- 
ties from LV distances, the SFI++ can also be exploited 
to study the RV and RL scaling relations of Sc galax- 
ies. In addition to well-characterized measures of L and 
V, a significant subset of SFI++ galaxies have homoge- 
neously measured as well as the radius r ^.R measured 
at the fii = 23.5 mag arcsec -2 isophote (jHavnes et all 
I1999U S07). Detailed studies using earlier compilations 
of the catalog explore and correct for internal extinctio n 
effects in these parameters (|Giovanelli et al.lll994 119951) . 



SFI++ disk galaxy scaling relations 
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Thus while the ideal approaches toward sample selection 
and parameter measurement may differ when compiling 
samp les for LV distances versus galaxy formation anal- 
yses (jPizagno et al.l [20071 : lAvila-Reese et al.l I2008L C07, 
see also £j5]), the size, quality and homogeneity of the 
SFI++ are unrivalled for both applications. 

In this paper, we present the LV, RV and RL scaling 
relations, their residuals and their scatter for a subset 
of the SFlH — h template and nontemplate samples. Since 
the LV relation has already been extensively studied by 
M06, we focus here on the RV and RL relations. We 
argue that 7-23.5 is superior to r^ as a measure of R, and 
present RL relations with observed scatters that are fac- 
tors of ~ 2.5—4 smaller than previously found. We derive 
detailed error budgets for all relations, and estimate the 
contribution of measurement errors to the observed scat- 
ters. We then discuss the applications RV and RL re- 
lations for measuring redshift-independent distances and 
constraining galaxy formation models, respectively. For 
clarity, we present the mathematical details related to 
the derivation of the SFI++ scaling relation parameters, 
their uncertainties, and the scaling relation error bud- 
gets in a pair of appendices. In all sections except £0 we 
adopt a value of Hq = IOO/1 -1 km s Mpc -1 for distance- 
dependent quantities. 

2. DATA 

We select our sample from the SFI++ (sec ££XJ) . Specif- 
ically, we include all SFI++ galaxies for which radii 
at the hi = 23.5 mag arcsec -2 isophote and disk scale- 
lengths have been homogeneously measured from /-band 
photometry (jHavnes et all Il999bl S07). We treat the 
template and nontemplate galaxies separately because 
(1) as discussed in £j4l peculiar velocities and incom- 
pleteness biases affect them differently and (2) differ- 
ent subsamples are appropriate for different applications. 
Throughout this paper, we refer to the 664/807 SFI++ 
template galaxies and 3655/4054 SFI++ nontemplate 
galaxies that meet our selection criteria as the template 
subsample and nontemplate subsample, respectively, and 
perform our analysis separately on each. 

2.1. Measurements, Corrections and Error Estimates 

Absolute /-band magnitudes and homogenized velocity 
widths for the SFI++ are presented in S07. The deriva- 
tion of these quantities, which draws on the work of the 
Cornell group ove r the past 15 years (IGiovanelli et al 



gnu 



2.1.1. Luminosities 



lM l [T995L fl997bt lHavnes et al.l I1999H ICatinella et al 



200a ISoringob et all 120051: ICatinella et al . 1 120071 M06. 
S07), is also summarized in that paper. Because these 
derivations are paramount to a rigorous computation of 
the scatter in the scaling relations presented here, we 
compile the relevant measurement, correc tion and un- 
certainty estimate equations in Appendices IA. II and IA.21 
We provide a brief description of these parameters and 
the adopted corrections in £12.1.11 and £|2.1.2I 

A number of disk size measures are also available 
for the template and nontemplate subsamples. In 
£ 12.1.31 we justify our choice of the radius at the /i/ = 
23.5 mag arcsec" 2 isophote as the disk size and give an 
overview of its derivation. The related mathematical de- 
tails are presented in Appendix I A. 31 

We follow the same prescription for computing dis- 
tances to SFI++ galaxies as in S07, which we describe 



Luminosities in the SFlH — h are expressed in terms of 
absolute /-band magnitudes, Mi. Their derivation is ex- 
plained in S07: apparent magnitudes extrapolated to 8 
disk scale-lengths are extracted from /-band photometry 
(jHavnes et al.lll999b[ ). They are corrected for Galactic 
extinction using the values of iSchleeel et alJ (Il998f) and 
internal exti nction using the relations of lGiovanelli et al. 
1994 , Il995l ). The type- dependent k-correction of IHan 
1992D is also applied. Absolute magnitudes are then 



computed using measured galaxy or cluster redshifts as 
summarized in £12.1.41 Uncertainties on Mj are com- 
puted by propagating (uncorrelated) measurement errors 
on the apparent magnitude and Galactic and internal ex- 
tinction corrections. The mathematical details of these 
computations, firs t presented in a series of papers by 
the Co rnell group (|Giovanelli et al.|[l997bl : lHavnes et alJ 
Il999bl M06, S07), are compiled in AppendixTATT" 

2.1.2. Rotation Velocities 

Rotation velocities in the SFI++ are expressed in 
terms of the logarithm of the velocity width, logVK, 
and are presented in S07. As explained in that paper, 
log W is derived from either single -dish HI profiles or op - 
tic al rotation curves (PR Cs) as in ISoringob et "ail §005) 
or ICatinella et al.l ([2005D . respectively. HI widths are 
corrected for instrumental broadening and turbulence, 
while ORC widths are measured from a parametric fit 
to the folded rotation curve and homog enized with the 
HI wi dths using the relations derived by Catinella et al. 
(|2007l) . All velocity widths are corrected for cosmologi- 
cal stretchi ng, and for inclination using measured /-band 
cllipticities (jHavnes et alJl!999bb IGiovanelli et al.lll997bl 
S07). Uncertainties on log 11^ are computed by propa- 
gating (uncorrelated) measurement errors on the veloc- 
ity width and cllipticity, as well as on the intrinsic ax- 
ial ratio of the disk. The mathematical details of these 
computations, firs t presented in a series of papers by 
the Co rnell group (IGiovanelli et al.|[l997U lHavnes et al] 
Il999bt ISpringob et all 12003 ICatinella et alJ 120051 120071 
M06, S07), are compiled in Appendix lA~2l 

2.1.3. Sizes 

As discussed in £JTJ deprojected disk scale-lengths rj are 
difficult to reliably measure despite their widespread use 
to construct scaling relations: the value of extracted 
from the surface brightness profile depends on the sub- 
jective process of choosing the profile's exponential re- 
gion ("marking the disk"), and extinctio n corrections ap- 
plied t o rd a re large and uncertain (e.g. IGiovanelli et all 
11994 119951 ). Given these effects, it is clear that 
the actual uncertainty in extinction-corrected is 
much larger than the ~ 15 % measurement errors that 
are generally assigned (e.g. iSchombert fc BothunI [19871: 
i Bvun fc Freeman! fl995l : Ide Jongl 119961 : IMacArthur et ail 
l2003HFathi et al.ll201oj i. 

We therefore adopt the isophotal radius 7-23.5 measured 
at the Hi = 23.5 mag arcsec" 2 isophote as the SFI++ 
measure of disk size. This quantity is straightforward 
to measure h omogeneously from the high-quali ty SFlH — |- 
photometry (jHavnes et all Il999bl see also iCourteaul 
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Fig. 1 . — Ratio of the measured dis k scale -length to the isopho- 
tal radius r'23.5 (computed using eq. IA23I I as a function of Mj for 
the nontemplate subsample. The points show the median ratio 
in equally populated absolute magnitude bins, and their errorbars 
represent the uncertainties on the position of the median deter- 
mined by bootstrapping. The solid line shows the best linear fit to 
the points. 

Il996| ). and does not require marking the disk. It is mea- 
sured at a location far enough from the bulge that the lat- 
ter does not contribute significant ly to the light. It s incli- 
natio n dependence was studied bv lGiovanelli et all (|1994 
|1995[) : they found that galaxies with Mj > — 21 mag are 
completely transparent at fii = 23.5 mag arcsec -2 , mak- 
ing the extinction correction also straightforward. 

We define the isophotal radius that we adopt for the 
SFI++ disk size as -R23.5 in units of /i -1 kpc, and ex- 
press it logarithmically. Wc presen t a detailed deriva- 
tion of logi?23.5 in Appendix IA.3I Briefly, the radius 
corresponding to the [ij = 23.5 mag arcsec -2 isophotc 
is measured f rom /- band photometry as described by 
lHavnes et afl(|1999bD . The same Galactic extinction and 
k-corrections adopted to compute Mj ( ^2.1.10 are then 
applied, and the resulting values are corrected for cos- 
mological surface brightness dimming and cosmological 
stretching. We then correc t for internal ex t inctio n us- 
ing the relations derived bv lGiovanelli et ail (|1995| ). and 
convert from angular to physical sizes using the measured 
galaxy or cluster redshifts a s explained in £12.1.41 

Following the approach of lGiovanelli et al.l (|1997bj) , we 
derive uncertainties for each value of logi?23.5 by propa- 
gating uncorrelated measurement uncertainties. Specifi- 
cally, we include uncertainties on the raw isophotal radii, 
disk scale-lengths, Galactic and internal extinction cor- 
rection, as well as the measured uncertainties on the disk 
cllipticity. 

Fig. [ T] sho ws the ratio of r<2 to ^3.5 (the latter given 
by eq. IA23[) as a function of Mj for the nontem plate 
subsample. As explained bv lGiovanelli et all (|1995l) . the 
dependence of rd/r 2 3.5 on Mr stems from the system- 
atic change of disk central surface brightnesses with this 
quantity. The solid line in Fig. [T| is the best linear fit to 
the median values indicated by the points: 

(-^-) = 0.672 + 0.0206(Af/-51og/i) . (1) 

We multiply -R23.5 by this relation to convert our disk 
sizes to scale-length units when comparing with previous 
studies (£J3J) and galaxy formation models f £15.21) . 

2.1.4. Distances 



Following M06 and S07, wc adopt different values of 
the CMB-framc rcdshift cz in the computation of Mr and 
log -R23.5 for the template and nontemplate subsamplcs. 
For the template subsample, "in" galaxies are assigned cz 
of their parent cluster measured by M06, while distances 
to "in+" galaxies are computed from cz measured for 
each galaxy assuming pure Hubble flow (see M06 for the 
definition of "in" and "in+"). For the nontemplate sub- 
sample, distances are computed from cz for each galaxy, 
assuming pure Hubble flow. 

Note that S07 compute peculiar velocities for nontem- 
plate SFI+4- galaxies, which they take as the offset be- 
tween Mr and the template LV relation from M06. How- 
ever, only part of the scatter in this relation stems from 
peculiar velocities (see £j4j, which implies that the values 
computed by S07 are upper limits. We therefore assume 
pure Hubble flow for our nontemplate subsample galax- 
ies, and treat peculiar velocities as a source of intrinsic 
scatter in the relations. 

3. SCALING RELATIONS IN THE SFI++ 

We now construct LV, RV and RL scaling relations 
for SFlH — h template and nontemplate subsamplcs. 

3.1. The LV , RV , and RL Relations 

Figures [2] and [3] show the LV, RV, and RL relations 
for the template and nontemplate subsamplcs defined in 
£J21 respectively. All scaling relations were modeled us- 
ing an orthogonal linear fitting method which takes into 
account measurements on x and y simultaneously. A re- 
liable estimate of the scatter e b s in the y-axis variable 
about the best-fitting linear relation is computed by ap- 
plying Tukey's biweight to derive a robust measurement 
of the dispersion^]. The best-fitting orthogonal relations 
and 2e b s intervals are plotted as solid and dotted lines 
in each panel of the figures, and the corresponding fit 
parameters are given in Table [T] 

The morphological corrections in cqs. [21 [3] and [4] have 
been applied to all of the relations in Figs. [21 and [3] (sec 
£13. 2j) . They are therefore representative of Sc galaxies 
specifically. We present the best-fitting linear relations to 
both the uncorrected and corrected relations in Table [TJ 
As a simple measure of the tightness of the relations, the 
Pearson correlation coefficient r is computed in each case 
and is also given in Table [lj as well as in the top-right 
corner of each panel in Figs. [21 and [3] 

The top cross in the bottom-right corner of each panel 
in Figs. [2] and [3] shows the median measurement uncer- 
tainty e mes in the sample, and the cross below it shows 
the median ± 3 times the median absolute deviation of 
that distribution. We discuss the sources of scaling rela- 
tion scatter in detail in £j4j for now, we point out that, as 
found in many previous studies, e Q h s easily exceeds the 
median e mes . 

There is a tight linear correlation between L and 
V in the SFI++; we refer the reader to M06 for 
an exhaustive discussion of the SFI4-+ LV relation. 
For the entire SFI++ template sample, M06 derived 
an incompleteness-corrected, morphology-corrected tem- 
plate LV relation parametrized by Mi — 5 log h — 

3 In the case of an outlier-free, Gaussian distribution this ap- 
proach reduces exactly to the classical standard deviation, but in 
the presence of strong outliers it provides much more stable results. 
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Fig. 2. — Scaling relations for the SFIH — h template subsample: (a) the luminosity-velocity (LV) relation, (b) the size-velocity (RV) 
relation and (c) the size-luminosity (RL) relation. The morphological corrections of eqs. [4] [2| and [3] have been applied in panels (a), (ft) 
and (c) respectively. In each panel the best orthogonal linear fit is plotted as a solid line, and the 2<r scatter 2e D (, s in the y-axis variable 
about that fit is delimited by dotted lines. The Pearson correlation coefficient of each relation is given in the top-left corner of each panel. 
The crosses in the bottom-right corner of each panel show the median measurement uncertainty in the sample (top) and the the median 
uncertainty ± 3 times the median absolute deviation of this quantit y. In (a), the dashed l ine and shaded region shows the SFIH — h template 
LV relation derived by M06 and the dot-dashed line the relation of Sakai & et al. (2000). A color version of this figure is available in the 
electronic edition of the Journal. 



-20.85 -7.85(logVF- 2.5). This relation and its median 
scatter are shown in Fig.[^K as the dashed line and shaded 
region. The M06 relation is slightly shallower than the 
one derived here due to small differences in sample com- 
position: our template subsample contains only ~ 80% 
of the galaxies used by M06 (sec f}2]). 

We note that M06 correct the SFI++ LV template for 
incompleteness, an effect that biases its slope low and its 
zero point high because we are more likely to detect the 
brighter galaxies near the flux limit of a survey (thus at 
low V). We do not correct the template subsample used 
here for incompleteness. As a cursory investigation of 
the impact of incompleteness on our sample, we fit the 
LV relation only to galaxies with logVF > 2.4, which 
M06 demonstrate are unaffected by incompleteness in 



the SFI++ template (see also ^4.1.2p . As expected, the 
re-fitted relation has a steeper slope and brighter zero 
point, pulling it farther away from the M06 result shown 
in Fig. [2k- Nonetheless, the slope and zero point of this 
re-fitted LV relation are consistent with those in Table [1] 
within their errors. At least for the LV relation, in- 
completeness therefore does not produce a strong bias 
relative to the uncertainties in our linear fits, in part be- 
cause of the smaller sample adopted here compared to 
M06 (and thus the larger statistical uncertainties). We 
discuss the effect of incompleteness on the scaling rela- 
tion scatter in SjU 

Taking similarities in sample selection, photometric 
band and parameter measurement into account, it is 
most straightforward to compare our nontemplate sub- 
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Fig. 3. — Same as Fig. [2] but for the SFH — h nontemplate subsample. The long-dashed lines and light (blue) shaded regions show the 
best-fitting relations and (2<r) scatter determined by C07, and the dash-dotted lines and dark (yellow) shaded regions show the best-fitting 
relations and scatter determined by P05. The parameters in these studies arc converted to our notation as described in the text. A color 
version of this figure is available in the electronic edition of the Journal. 



sample relations with the results of P05 and C07. The 
best-fitting relations and 2a scatters from those studies 
arc ovcrplotted on our nontemplate relations in Fig. [3l 
To carry out the comparison, luminosities are converted 
to magnitudes assuming Mqj = 4.19, and we correct for 
the slight difference in photometric band of the P05 anal- 
ysis using log Li = log Lj + 0.036 (C07). We change from 
rotation velocities to velocity widths by multiplying by 
two. Finally, we convert from the scale-length units used 
in P05 and C07 to the isophotal radius units adopted here 
using the relation in eq. [TJ Considering the differences in 
sample selection and parameter estimation, there is rea- 
sonable agreement between the best-fitting LV relation 
and scatter found here and those of P05 and C07 (dash- 
dotted and long-dashed lines in Fig. El respectively). 

We also find strong correlations between R and both V 
and L in the SFI++. Fig. [3] shows that when converted 
to our units, the best-fitting C07 relations agree very well 



with ours. This is expected given the significant overlap 
between the samples, but reassuring to confirm because 
of the different parameters and corrections adopted in 
both studies. By extension from the discussion in C07 
and more recent studies, there is broad general agree- 
ment between our best-fitting relations and others in 
the literature (Ide Jong fc Lacevl 120001: iShen et all 120031: 
iGraham fc Worleyll2008l : lFathi et al.ll2010l) . As suggested 
bv lDutton et al.l (|2007f ). the systematically larger scale- 
lengths found by P05 may stem from differences in galaxy 
selection and parameter derivations relative to other 
studies, and is likely exacerbated by their small sample. 

In contrast to the similarity between our best-fitting 
RV and RL slopes and zero points and those in the lit- 
erature, the high correlation coefficients of our relations 
(Table [1]) provide the first indication that they are sig- 
nificantly tighter than has been previously reported: for 
example, we measure r = 0.84 and r = 0.97 for the RV 
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and LV relations in the nontcmplatc subsample, while 
C07 measure an average of r ~ 0.65 for the analogous 
relations. 

Fig. [3] illustrates how the scatter in our best-fitting RV 
and RL relations (particularly the latter) is significantly 
smaller than that reported by P05 and C07. Focussing on 
the more widely studied RL relation, we find an average 
scatter of e & s = 0.05 log(/i _1 kpc) for both the template 
and nontemplate subsamples (Table[l}. Converted to our 
units, C07 report a scatter of <Ji og R ~ 0.14 log(/i _1 kpc) 
using a sample that is similar to and overlaps with the 
SFlH — h in many respects. P05 emphasize their small, 
well-characterized measurement uncertainties in their ho- 
mogeneous sample of 81 galaxi es, and report an RL scat - 
ter identical to that of C07. lAvila-Reese et all (|2008l ) 
state that their heterogeneous sample of 76 galaxies 
spans a broad range of morphological types and sur- 
face brightnesses; they find cji „ „r ~ .20 log(fe~ 1 kpc) for 
the RL relation. IShen et all |2003[) study the RL re- 
lation for a statistically complete sample of ~ 100 000 
late-type galaxies from the SDSS, selected using cuts 
in light concentration and Sersic index. They report 
o~iogR. ~ 0.13 log(/i _1 kpc) at the high-luminosity end of 
their relation (late and e arly types are likely mixed at the 
low- luminosity end; see IGraham fc Worlevl I2008D . We 
therefore find that relative to studies adopting a wide 
range of sample selection philosophies, the average scat- 
ter in the SFlH — \- RL relation is factors of 2.5 - 4 smaller 
than previously reported. 

We attribute the significantly smaller scatter in our RV 
and RL relations to our use of homogeneously measured, 
extinction-corrected isophotal radii as disk sizes instead 
of the scale-lengths used in the above studictQ. As ex- 
plained in £12.1.31 there are good reasons to expect that 
i?23.5 can be more reliably measured and corrected than 
r<j, and therefore that it should produce tighter scaling 
relations than the latter. We carry out a detailed analysis 
of the scatter in the RL relation in fj4] 

Fig. |4] shows the nontemplate RL relation constructed 
in exactly the same manner as in Fig[3t, except that we 
use the SFI++ disk scale- length r<i, co rrected for inter- 
nal ex tinction using the prescription of IGiovanell i et all 
(|1995f ). instead of i?23.5- The values of r and e bs of that 
relation (including all galaxy types) are given in the top- 
left corner of the plot. This version of the RL relation 
has a similar r to that found by C07, and e b s that is 
comparable to the scatters reported in the studies de- 
scribed above. Fig. [4] confirms: (1) the speculation by 
C07 that RV and RL relations constructed using are 
less robust than the LV relation because is not re- 
liably measured, (2) that extinction-corrected 7*23.5 ar e 
superior to extinction-corrected r^ when building scaling 
relations in that the former parameter yields significantly 
stronger correlations, and (3) that the small scatter in 
our RL relations is not caused by s elective "pruning" 
of hig h-scatter points in the SFI++ (jAvila-Reese et all 
2008). 



4 IShen et al.l 1 120031) use Petrosian and Sersic half-light radii in 
their study, not r&. Tests using the radius encompassing 83% of the 
/-band light for SFlH — h galaxies show that integral measures also 
produce scaling relati ons with more scatter than presented here 
(Saintongc et al. 2008), presumably because of the inclusion of the 
bulge light and lack of extinction correction. 
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Fig. 4. — The RL relation for the nontemplate subsample, using 
the extinction-corrected disk scale-length (r^) as the size indicator 
instead of -R23.5- The data are broken down by morphological type, 
with lines and symbols described in Fig. \E\ The best orthogonal 
linear fit to the entire sample is given by the thick solid line, and the 
scatter 2e (, s about that fit is shown by the thick dotted lines. The 
values of e a i, a and the Pearson correlation coefficient r computed 
for the entire subsample are in the top-left corner. A color version 
of this figure is available in the electronic edition of the Journal. 

3.2. Morphological Dependence 

As discussed by M06, the SFI++ contains galaxies 
with a broader range of morphological types than ear- 
lier incarnations of the catalog, affording a more thor- 
ough investigation of scaling relations as a function mor- 
phological type. Nonetheless, the histograms in Fig. \5]p 
illustrate that the SFI++ is dominated by Sc galax- 
ies. As such, it is difficult to disentangle sample biases 
from physical effects in any scaling relation differences 
that are found. We t herefore follow the approach of 
IGiovanelli et al.l (|1997bf ) and M06, and adjust the best- 
fitting scaling relations determined for galaxy types Sa 
and Sb to match that of the best-fitting Sc+Sd relations. 
We use the larger nontemplate subsample for this exer- 
cise, and apply the derived correction to both the tem- 
plate and nontemplate subsamples. 

The computed slopes, zero-points and scatters of the 
scaling relations to which these morphological corrections 
have and have not been applied are given in Table [TJ 
There is little difference between the properties of the 
uncorrected and corrected relations: this stems from the 
high Sc fraction in the SFI++. For this reason, none 
of the conclusions in our study are affected by our ap- 
plication of the correction. The morphological depen- 
dence of the LV relat ion has been extensively discussed 
in th e literature (e.g. iRobertsI [19781: lAaronson fc Moulc 



TM a IRubin et all 11985): IGiraudl Il986at IPierce fc TuTiv 



1988t iKannappan et al.N2002L C07). That of the SFI++ 



LV relation in particula r has already been addressed by 
IGiovanelli et al.1 ()1997bl ) and M06. We therefore focus 
on the morphological dependence of the RL relation, and 
provide corrections for the LV and RV relations using 
the same method. 

Fig. [5jt shows the RL relations for nontemplate galax- 
ies of different morphological types, as classified in Ta- 
bles 2 and 4 of S07. Note that for clarity, we only plot 
a random subset of 100 galaxies of each morphological 
class. The lines in Fig. [5jt show the best orthogonal lin- 
ear fits to galaxies of each morphological type. We detect 
a clear trend: at a given L, early-type spirals have smaller 
R than late-type spirals, and the slope of the early-type 
RL relation is steeper. This trend is qualitatively simi- 
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TABLE 1 

Orthogonal Fit Parameters for SFH 



Scaling Relations 



Relation 




Template 


Non-template 








a ±da 


b ±db a r a±da 


b±bb 


cr 


r 


NO MORPHOLOGICAL CORRECTION 


LV a : 
RV b : 
RL C : 


-20.83 ±0.17 
0.911 ±0.037 
0.763 ± 0.028 


-8.79 ±0.04 0.55 0.91 -20.77 ± 0.07 
1.281 ±0.011 0.11 0.84 0.910 ±0.015 
-0.159 ±0.001 0.05 0.96 0.763 ±0.016 


-9.07 ±0.02 
1.357 ±0.004 
-0.165 ±0.001 


0.54 
0.11 
0.05 


0.87 
0.83 
0.97 


WITH MORPHOLOGICAL CORRECTION 


RV h : 
RL C : 


-20.88 ±0.15 
0.927 ±0.038 
0.774 ±0.029 


-9.24 ±0.04 0.55 0.91 -20.81 ± 0.06 
1.353 ±0.011 0.11 0.85 0.922 ± 0.019 
-0.160 ±0.001 0.05 0.97 0.773 ± 0.012 


-9.42 ±0.01 
1.407 ± 0.005 
-0.164 ±0.000 


0.53 
0.11 
0.05 


0.88 
0.84 
0.97 



a The LV relation is parametrized as Mj — 5 log h = i 
b The RV relation is parametrized as logR [h —1 kpc] 
c The RL relation is parametrized as log R [h^ 1 kpc] = 



■b(\ogW -2.5). 
i + b(logW - 2.5). 
, + 6(Afr+20). 



lar to that found bv lShen et al.l (|2003l ) for their complete 
sample. Based on the orthogonal fits in Fig. [5k, we apply 
the following additive factor to L in the RL relation so 
that the mean relation for early types matches that of 
the later types: 



A 



RL 




0.024(M/ + 20) for types SOa/Sa 
0.008(Mj + 20) for types Sab/Sb 
for later types. 



(2) 

Note that this correction is only applied to L in the RL 
relation; L for individual Sa and Sb galaxies in the LV 
and RL relations of Figs. [2] and [3] are thus slightly dif- 
ferent by virtue of the morphological corrections applied. 
We include a 15% uncertainty in this correction in the 
RL relation error budget (see Appendix iBj). 

We find a lesser dependence of the RL relation param- 
eters on morphology than reported by C07. This differ- 
ence stems from our definition of R: the lines in Fig. U 
show the best-fitting relations constructed with for dif- 
ferent morphological types, where we recover a strong de- 
pendence. The systematic variation of (rrf/r23.s) in Fig.[T] 
can explain at least part of the difference in morpholog- 
ical dependence between these two indicators: relative 
to rd, r 2 3.5 is systematically larger for brighter galaxies, 
which are more likely to be Sa's according to Fig. [SJd. 
This works to offset the strong trend found when is 
used to construct the RL relation. However, we can- 
not rule out systematic effects in how is measured in 
Sa's relative to Sc's, such as a bias in the manner that 
Sa disks are marked given their more substantial bulges 
(see g2Tp. 

We derive the morphological corrections to the RV and 
LV relations in an analogous manner to that for the RL 
relation: 

f 0.77 - 0.25(log W - 2.5) for types SOa/Sa 
A RV = I -0.13 + 0.06(logPF - 2.5) for types Sab/Sb 
[ 0.0 for later types. 

(3) 

f -1.24 + 0.31(log W-2.5) for types SOa/Sa 
A LV = I 2.54- 1.01(logVF- 2.5) for types Sab/Sb 
{ 0.0 for later types. 

(4) 

The small difference between the correction in eq. [4] and 
that from M06 is likely due to the fitting techniques (bi- 




-20 -21 
M, - 5 log h 

Fig. 5. — a) RL relation for the nontemplate subsample, subdi- 
vided by morphological type. The lines represent the best orthog- 
onal linear fits to the different subsets. To avoid overcrowding, 
only 100 randomly drawn galaxies of each morphological type are 
plotted. The definition of each morphological class based on the 
T-types of S07 is given in the upper-left corner, b) Normalized his- 
tograms showing the luminosity distribution of each morphological 
type, with the total number of galaxies of each type given. A color 
version of this figure is available in the electronic edition of the 
Journal. 

sector fits in M06 and orthogonal fits with x and y errors 
here) and to sample selection: we calculate our correc- 
tions using the larger nontemplate subsample, while M06 
use their incompleteness-corrected SFIH — h template rela- 
tion. 

3.3. Scaling Relation Residuals 

The residuals of the scaling relations provide useful 
checks on data quality as well as important insight into 
the nature of the the relations themselves. Here, we ex- 
ploit the large size of the nontemplate subsample to ex- 
amine the residuals of the scaling relations presented in 

Fig. [S] plots the data-model residuals of the LV, RV 
and RL relations as a func tion of J-band centra l surface 
brightness fig (derived as in lHavnes et aLlll999bf ) and the 
third scaling relation parameter. In all panels, the resid- 
ual is computed in the variable on the y-axis in Fig. |3l 
We find no significant correlations between any of the 
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Fig. 6. — Data - model residuals of the scaling relations for the 
nontcmplatc subsamplc, plotted as a function of central /-band 
surface brightness (a, b, c) and the third scaling parameter (d, e, 
f), namely R for the LV relation, L for the RV relation and V for 
the RL relation. In all panels, contours show the distribution of 
all datatpoints. The red points show the median value in equally 
populated bins, with errorbars representing the la spread in the 
distribution of values. 

scaling relation residuals and these parameters. We have 
also checked for correlations with variables such as incli- 
nation that might imply biases in our samples, and find 
none. 

The lack of correlation between the residuals of the 
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ally i nterpreted as evidence in favor of sub-maxim a l disks 
(e.g. iBershadv etafl 120101 ). but iDutton et all (|2007l ) 
demonstrate that a variety of factors can alter the sur- 
face brightness dependence of the LV relation (se e also 
iFirmani fe Avila-Reesdl2000l: Ivan den Bosch! [2000) . 

On the other hand, the lack of correlation between the 
residuals of the RL and RV relations and fio is not typ- 
ical of past survey results. The reason is again our def- 
inition of R: Fig. [7] shows that we recover a strong sur- 
face brightness correlation when the extinction-corrected 
scale-length is used as the measure of disk size in the 
RL relation, with a slope consistent with that reported 
by C07. In analogy to the morphological dependence 
discussed in £13. 2| the change in (r^/ V23.5) as a func- 
tion of L in Fig. Q] can also explain the lack of cor- 
relation when -R23.5 is used to construct the RL rela- 
tion: brighter, higher surface brightness disks have lower 
which reconciles Fig.[7]and Fig. [BJ;. However, 
it is also possible that systematic effects related to the 
measurement of and /io (the latter an extrapolation of 
the exponential disk defined by the former) also influence 
the trend in Fig. [7] 

Finally, we investigate the correlation between the scal- 
ing relation residuals, shown in Fig. [8] The axis labels in 
Fig. [8] reflect the definition of the residuals: for example, 




Fig. 7. — Residuals of the RL relation for the nontcmplatc sample 
constructed using the extinction corrected disk scale-length instead 
of ^?23.5 (Fig. [4}, plotted against the /-band central surface bright- 
ness. The best-fitting linear relation is shown by the solid line, and 
is given in the bottom-right corner. 

ALR on the y-axis of Fig. [5^ are the residuals of the 
RL relation, computed as data - model residuals in L, 
whereas ARL on the y-axis of Fig. [8]o are the residuals 
of the RL relation, computed as data - model residuals 
in R. FigsEKjIHb and[8fc therefore correspond to the LV, 
RV and RL relation in differential form. 

As in P05 and C07, we find little evidence for a correla- 
tion between the LV an d RL relation residuals in Fig.[5jD. 
iCourteau fcRixl (fl999l) argue that this lack of correlation 
suggests that even high surface brightness galaxies have 
subma xima l disks. Howev e r, the models of lGnedin et al.1 
(|2007f) and IDutton et all (|2007t) demonstrate that fac- 
tors such as bulge formation, stellar and gas fractions, 
scatter in halo properties and disk mass-to-light ratios 
also influence the behavior of the residuals. We find a 
weak positive correlation between the RL and RV rela- 
tion residuals, and a strong positive correlation between 
the LV and RV relation residuals, whose slopes are in 
rough agreement with those of the actual LV and RV 
relations, respectively. As explained in C07, these re- 
sults are qualitatively consistent with little scatter in L 
and uncorrelated scatter in V and R. The error budget 
analysis in ^supports these conclusions for the measure- 
ment errors in V , L and R, but the substantial intrinsic 
scatter in all of the relations requires detailed modeling 
in order to fully characterize the origin of the residual 
correlations. 

4. ERROR BUDGETS OF THE SCALING RELATIONS 

An understanding of the observed scaling rela- 
tion scatter e f, s is essential for reconstructing large- 
scale structure from extrac t ed peculiar vel o cities (e.g. 
Sandage fc Tammannl 119751 : IGiraudl Il986bt IDale et al.l 



1999aD as well as for modeling the origin of the rela- 



tions themselves (e.g. IDutton et al.l 120071 : iGnedin et al.l 
120071 see jfl]). In the same spir it as the analyses carried 
out bv lGiovanelli et al.l ()1997bD and M06 for the SFI++ 
template LV relation, we compute the error budgets for 
both the template and nontemplate LV, RV and RL 
relations presented in fJ3] 

The total observed scatter in the scaling relations, e b s , 
is a combination of the measurement uncertainties e mes 
and ei n t, the intrinsic scatter: 

2 _ 2 I 2 fr~\ 

For each relation, e mes is computed from the measure- 
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Fig. 8. — Correlations between the residuals of the LV, RV and RL relations for the SFH — (- nontemplate subsamplc. The axis labels 
denote how the residual is defined in each case: for example, A.LR in (a) are data - model residuals in L computed from the RL relation, 
whereas ARL in (ft) are data - model residuals in _R computed from the RL relation. Panels (a), (6) and (c) therefore represent the LV, 
RV and RL relations in differential form, respectively. In each panel, the Pearson correlation coefficient of the residuals is in the top-left 
corner, the solid lines are the best-fitting linear relations to the residuals, and their slope is in the bottom-right corner. The dashed lines 
represent the slope of the corresponding scaling relation from Fig. \3\ 
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Fig. 9. — Error budget of the LV relation for the SFI++ template 
(a) and nontemplate (b) subsamples. Data points are the median 
observed scatter (e f, s ) as a function of V, and the error bars are 
the uncertainty on the position of the median. The separate con- 
tributions to the total scatter are: the measurement error on L 
(e^, dotted line), the measurement error on V (\b\ e\y, dash-dotted 
line), which add in quadrature to produce the light solid line (i.e. 
neglecting the contribution of corre lated errors). The total mea- 
surement error e mea = ej,y (eq. IB2I lower dark solid line) accounts 
for correlated measurement uncertainties and, for the nontemplate 
subsample (ft), the contribution of a peculiar velocity of amplitude 
Vpec = 300 kms - 1 (eiv.Vpec, dashed line). The solid line that 
goes through the points at logW > 2.4 is the sum in quadrature 
of Ely and of a constant intrinsic scatter ei nt . A color version of 
this figure is available in the electronic version of the Journal. 

merit uncertainty on the two scaling parameters and their 
covariance. We also include a contribution from the mor- 
phological corrections of £13.21 The mathematical details 
of the derivation of e mes for the three scaling relations 
can be found in Appendix [B] 

4.1. The LV and RV relations 

Figures and [TU] show the total scatter of the LV 
and RV relations as a function of V for both SFI++ 
subsamples. We note that our estimates of e mes (dash- 
dotted and dotted lines) for each varia ble are simi- 
lar to that reported in o ther studies (e.g. IDutton et al.l 
l2007t iGnedin et al.l 120071) , and that accounting for cor- 
related uncertainties does not significa ntly change e mes 
(lower dark and light solid lines; see also lGiovanelli et al.l 
Il997bl) . 



Fig. 10. — Error budget of the RV relation for the SFI++ tem- 
plate (a) and nontemplate (ft) subsamples. Lines and symbols are 
as in Fig. [9] A color version of this figure is available in the elec- 
tronic version of the Journal. 

In all cases, the amount of scatter is almost constant 
at loglF > 2.4, but increases at lower V. Above this 
threshold, e b s can be very well reproduced by adding 
a constant intrinsic scatter cm in quadrature with e mes 
(elv and €rv] the red line in each panel). 

The value of Ei n t required to match e (, s at logVF > 
2.4 for the template and nontemplate subsamples are 
broadly consistent, and we attribute the slightly larger 
value in the template case to a small peculiar velocity 
effect in that subsamplc (sec I4.1.1j) . Using the values 
derived from the nontemplate subsample, which bene- 
fits from much better statistics, we adopt intrinsic scat- 
ters eLV.int = 0.35 ± 0.01 mag and £Rv,int = 0.084 ± 
0.001 log(h _1 kpc) for the LV and RV relations, respec- 
tively. 

Our value of thv,int is similar to that invoked by M06 
in their error budget of the SFI++ template LV rela- 
tion, although they fit for a width-dependent scatter in- 
cluding points with logll 7 < 2.4. That solution is per- 
fectly consistent with our error budget of Fig. but 
the cleaner view provided by the large nontemplate sub- 
samplc (Fig. [Sb) makes us favor the constant scatter ap- 
proach. 

As for the deviation from constant intrinsic scatter 
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at low V in Figs. H] and QJJ1 we identify four possible 
causes: peculiar velocities, sample incompleteness, in- 
creased measurement errors and increased intrinsic scat- 
ter. We investigate these possibilities in turn below. 

4.1.1. Peculiar velocities 

Peculiar velocities introduce scatter in both the LV 
and RV relations. Because distances to nontemplatc 
galaxies are derived assuming pure Hubble flow, the ef- 
fect will be strongest in that subsample. While peculiar 
velocities in clusters are larger than in the field, the tem- 
plate sample is mostly free of peculiar velocity-induced 
scatter, because we adopt cluster distances for all their 
members (see §2.1.4j) . 

The median redshift cz of the galaxies in the nontem- 
plate subsample is a strong function of V. Since a galaxy 
at redshift cz with peculiar velocity V pec will scatter away 
from the mean LV and RV relations in proportion to 
log(l + Vp ec /cz) (sec Appendix IB) , the effect of peculiar 
velocities on the scatter in these relations is strongest at 
low V , where cz is smallest. 

To assess the impact of peculiar velocities, we compute 
the scatter produced in the LV and RV relations by a 
characteristic V pec = 300 km s -1 and function of V 

in the nontemplate subsample. The result is illustrated 
by the dashed lines in Figs. [9Jd and ITOb . At the low- 
V end of the nontemplate relations, peculiar velocities 
can constitute the largest source of scatter. The lower 
dark solid line in these figures adds this contribution to 
the scatter in quadrature with e mes . The resulting is a 
steepening of the total measurement error curve at low V. 
However, it is clear that peculiar velocities alone cannot 
explain the increase in e b s in the LV and RV relations 
at low V. 

4.1.2. Sample incompleteness 

In addition to biasing scaling relation slopes and zero 
points (see EQ.ip . statistical incompleteness of a sample 
can also affect their scatter. Cluster samples are often 
used for scaling relation studies because they are less 
affected by incompl eteness than flux - limited field sam- 
ples. Nonetheless, iGiovanelli et al.l (|1997bl ) and M06 
have shown that the SFI++ template subsample suffers 
from incompleteness at logW < 2.5 (e.g. Fig. 4 in M06). 

The nontemplate subsample certainly also suffers from 
incompleteness. However, the union of several disparate 
samples to create the SFI++ makes it impossible to de- 
rive the selection function of the survey and quantify 
its incompleteness (see $TJ and the discussion in S07). 
It is nonetheless reasonable to assume that the qualita- 
tive impact o f incompleteness mirrors that determined by 
M06 (see also lde Jong fc Lacevll2000l ) . bmce we are more 
likely to observe larger/brighter galaxies near the survey 
limits, the scatter at the low-V end of the relations is 
reduced. However, we observe an increase in e t, s . Thus 
unless incompleteness in our SFI++ subsamples behaves 
in the opposite manner from that determined by M06 for 
the SFI++ template LV relation, it is unlikely to explain 
the increase in e b s at logVF < 2.4 in Figs. 1^1 and [TU1 

4.1.3. Increased scatter in V 

Contrary to the LV and RV relations, e b s in the RL 
relation does not increase at low L (see ^4.2l and Fig. [TTj) . 



For this reason, we associate the behavior of the scatter 
at low V in the LV and RV relations to V itself. Here, 
we investigate the possibility that there is an additional 
contribution to the scatter in V for which we have not 
accounted in our error budget. 

Part of the increased scatter could come from mea- 
surement uncertainties specific to galaxies with low V. 
Galaxies with narrower HI lines are more prone to 
"catastrophic" measurement errors due t o, for example, 
line a symmetries and turbulence effects (jSpringob et al.l 
12005V or uncertainties in inclination. While it is im- 
possible to account quantitatively for the effect of these 
catastrophic errors on our scaling relations, it is plausi- 
ble that they produce a total measurement error function 
{clv, £rv) that is steep at the \ow-V end. Nonetheless, 
it seems unlikely that the median fractional uncertainty 
on log W at log W ~ 2.2 is underestimated by a factor of 
4 relative to that computed in our budget, and thus that 
measurement errors alone explain the increase in e f, s . 

We are therefore left with the possibility that the scat- 
ter at low V stems from an astrophysical effect. For 
example, galaxies with low V have lower masses, and are 
therefore more susceptible to HI disk stripping processes. 
Sinc e these galaxies tend to have rising rotation curves 
(e.g. iCatinella et~aT]|2006t) , measuring V at different radii 
because of HI disks stripped to various degrees will in- 
troduce scatter. However, the slope of low-mass galaxy 
rotation curves is not nearly steep enough to explain the 
entirety of the increase in e b s . 

It is possible that the behavior of e b s in the LV and 
RV relations stems from an increase in the intrinsic scat- 
ter at low V related to the properties of their stellar disks 
or their parent dark matter halos. While this is a tan- 
talizing possibility at the outset, it is difficult to "tune" 
galaxy formation models of the type described in £ 15.21 to 
increase the scatter in the LV and RV r elations with- 
out in creasing that in the RL relation (e.g. lDutton et all 
2007). It is therefore not obvious that the increase in 
e bs at logPF < 2.4 in the LV and RV relations can be 
attributed to intrinsic scatter, although a full suite of 
galaxy formation models is needed to address this issue. 

4.2. The RL relation 

The average e b s for the RL relation is significantly 
lower than that of the RV relation (Table [TJ). This comes 
in part because the measurement errors on L introduce 
less scatter than those of V, as pointed out by C07 and 
clearly seen in the error budgets of Figures [10] and [TTJ 
However, our analysis shows that the relation also has 
significantly less intrinsic scatter: in regions at high V 
and high L where the scatter is approximately constant, 
we measure env,mt — 0.084 ± 0.001 log(/i _1 kpc) for the 
RV relation but e RL ,int = 0.034 ± 0.001 log(/i _1 kpc) for 
the RL relation. 

The overall behaviour of the scatter in the RL rela- 
tion is also different: contrary to the increases found for 
the LV and RV relations, we find a reduced scatter for 
Mi — 5 log h > —20 in the RL relation. We can rule out 
peculiar velocities affecting the scatter since they move 
points almost exactly along the RL relation (see Ap- 
pendix |B|) . Given the discussion in ^4.1.21 incompleteness 
is the most likely culprit. Since some of the data included 
in the SFIH — h stem from aperture diameter-limited sur- 
veys (especially among the fainter, more nearby galax- 
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e m =0.034±0.001 " 
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Fig. 11. — Error budget of the RL relation for the SFI++ tem- 
plate (a) and nontemplate (ft) subsamples. Lines and symbols are 
as in Fig. [9] In this case we do not include the contribution of 
peculiar velocities to the total scatter, because as discussed in Ap- 
pendix [Bj they move points almost exactly along the RL relation. 
A color version of this figure is available in the electronic version 
of the Journal. 

ies), we expect incompleteness to reduce the scatter at 
the faint end of the RL relation in analogy the LV rela- 
tion behaviour determined by M06 for the SFI++ tem- 
plate. 

5. APPLICATIONS OF THE RL AND RV RELATIONS 

In 8 j3.ll we present RV and RL relations that ex- 
hibit significantly smaller scatter than has been previ- 
ously reported, which we attribute to our adoption of an 
inclination-corrected isophotal radius as the size param- 
eter R in these relations rather than a disk scale-length. 
Here, we discuss the implications of this reduced scatter 
for two common scaling relation applications: redshift- 
independent distances and galaxy formation models. 

While it is obvious that the ideal samples to be used 
as distance indicators and galaxy formation constraints 
differ, it is equally obvious that these ideal samples are 
not yet available. The use of the scaling relations as ei- 
ther distance indicators or galaxy formation constraints 
thus requires finding an acceptable middle ground with 
respect to sample selection. It is in this sense that sub- 
sets of the SFI++ arc relevant in both the distance indi- 
cator and galaxy formation contexts. The SFI++ tem- 
plate sample was explicitly constructed to measure LV 
distances (M06), and we discuss in ^5.11 the potential of 
exploiting our template RV relation for the same pur- 
pose. The SFlH — h nontemplate galaxies span a broader 
range of properties, and were not selectively edited for 
cosmic flows or other purposes. The mix of morphologi- 
cal types in the SFI++, as with most other large catalogs 
of similar quality, is not representative of that in the lo- 
cal Universe. On the other hand, there is no reason to 
suspect that the SFI++ does not reflect the properties of 
Sc galaxies in that volume. With this caveat in mind, we 
explore in i)5.2l thc implications of the small scatter in the 
nontemplate RL relation for galaxy formation models. 

Note that in this section, we adopt the WMAP5 value 
of H = 71.9^2 ? km s" 1 Mpc^ 1 (IDunklev et all 12001 
for concretcness. 



5.1. The RV Relation as a Distance Indicator 

Since R depends on distance while V does not, the RV 
relation could in principle be used to produce distance- 
independent rcdshifts just like the widely-used LV rela- 
tion. This application has received little attention so far, 
because the RV relation typically has significant scatter. 
However, our RV relation is nearly as tight as the LV re- 
lation, and has ~ 1.5 times less scatter than found in pre- 
vious studies (see Figs. [2] and £j3j . Since disk sizes can be 
more straightforward to compute than magnitudes (and 
sometimes even more reliable), there is obvious interest 
in using the RV relation as a distance indicator. 

As a proof of principle, we use the sample of 17 SFI++ 
galaxies with a distance measurement from Cepheid vari- 
ables compiled by M06 to evaluate the distances derived 
from the SFI++ RV relation. We use the template sub- 
sample RV relation for this exercise, since it is less sus- 
ceptible to peculiar velocities and can in principle be 
corrected for incompleteness. The template sample in 
M06 was used to independently calibrate the zero point 
of the LV relation and to provide an estimate of Hq. 
Their value of Hq = 74 ± 2 (random) ± 6 (systematic) 
km s _1 Mpc -1 compares well with the WMAP5 result of 
H = 71.9±|? km s" 1 Mpc" 1 (jDunklev et al.|[2009h . The 
systematic component of the uncertainty in the M06 es- 
timate stems largely from the uncertainty of the Cepheid 
dista nce zero point, not the LV relation (iSakai et al.l 

12551 . 

Fig. HU shows the RV relation for the Cepheid sample, 
overlaid on the template RV relation normalized to the 
WMAP5 value of Ho- Since most of the galaxies in the 
Cepheid sample have large rotation velocities, it is not 
possible to reliably fit the slope of their RV relation. 
Instead, we adopt the slope of the template RV relation, 
but determine the zero point of the Cepheid sample RV 
relation independently. The ler region around the RV 
relation with this Cepheid-derived zero point is shown in 
Fig. [12] as the shaded region, which overlaps comfortably 
with that derived from the template subsamplc (solid 
line) . 

Reversing the problem, the value of Hq can be deter- 
mined by matching the zero points of the SFlH — \- tem- 
plate relation and that of the cosmology-independent 
value derived from the Cepheid sample. We find a value 
of Ho = 72 ± 7 km s _1 Mpc -1 , where the errors come 
in part from the intrinsic scatter of the RV relation but 
mostly from the uncertainty on the Cepheid distances. 
Comparing the M06 estimate of Ho with this value, we 
conclude that the LV and RV relations provide similar 
constraints on this parameter. 

Finally, the distances to the galaxies in the Cepheid 
sample are computed from the template subsample LV 
and RV relations in H3.1{ The relative differences be- 
tween the Cepheid- and scaling relation-derived distances 
are shown in Figure 1131 The typical errors on the LV- 
and i?V-derived distances are obtained by taking the me- 
dian absolute deviation of these relative differences, and 
are 14.8% and 24.7%, respectively. These errors are in 
line with the typical er rors of ~ 15—20% gen erally quoted 
for LV distances (e.g. lTullv fc Piercdl200"oT) . 

The above exercises suggest that even given the tight 
RV relation presented here, the LV relation produces 
more accurate redshift-independent distances. However, 
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Fig. 12. — RV relation for galaxies with distance measurements 
from Cepheid variables (points). The filled contours shows the 
template RV relation normalized to Ho = 71.9 km s _1 Mpc - x , 
and the solid line is the best fit to these data. The shaded region 
shows the la uncertainty in the zero point of the Cepheid RV 
relation when its slope is held fixed to the template RV relation 
value. 
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Fig. 13. — Relative difference between the distances from the 
Cepheid variable measurements and the distances derived from the 
SFIH — h scaling relations (the RV and LV relations, represented as 
filled circles and open squares, respectively). 

in situations where accurate luminosities are not avail- 
able, a well-defined RV relation can be substituted with 
only moderately larger uncertainties. We note that a 
thorough analysis of the RV relation as a distance in- 
dicator requires a correction for incompleteness effects. 
This is clearly feasible for the SFI++ template sample 
(jGiovanelli et al.lll997bl M06), but we defer this task to 
a future paper. 

5.2. The RL Relation and Galaxy Formation Models 

It is clear that a simultaneous application of galaxy 
formation models to the sca ling relations d erived here, 
like that ca r ried out by iDutton et al.l (|2007[ ) and 
iGnedin et al.l (|2007l ). is beyond the scope of this paper. 
Nonetheless, a consideration of the SFI++ RL relation 
scatter in the context of a simple model demonstrates 
the potential of our relations to place new constraints on 
galaxy formation. The RL relation is ideal in this context 
because peculiar velocities move points almost exactly 
along the relation, and therefore do not contribute to the 
scatter (see Appendix |B|) . We use the larger nontemplatc 
subsamplc for this exercise, with the morphological cor- 



rection derived in t]3.2l applied. This relation therefore 
constrains the formation and evolution of present-day Sc 
galaxies. 

5.2.1. Model Details 

We consider a IMo et all (|1998t hereafter MMW98) - 
style model of the SFI++ RL relation and its scatter, 
where a self-gravitating exponential disk is embedded in 
an NFW halo. Specifically, we eliminate the halo velocity 
V200 from their eqs. 16 and 28 to produce the predicted 
relationship between the disk luminosity Ld and scale- 
length Rd of the stellar mass distribution: 



Rd — 



G 



200V2F n 2 



1/3 T l/3 



A' 



1/3 



LTfc 1/2 f- 



(6) 



In eq. El is the dis k mass-to-ligh t ratio . We adopt 
the relation derived by IDutton et al.l (l2007f) us ing the 
population synthesis models of I Bell et al.l ( 20031 ) and a 
"diet Salpeter" IMF: 



log 



T d ,i 



M o /L ,i 



= 0.172 + 0.144 log 



10 



10.3 



'0,1 



(7) 

The parameter vrid in eq. |B]is the ratio of the disk mass 
Md to the halo mass M2oo- Few constraints exist on 
or on its scatter at fixed M200, b ut most spiral galax y 
mass models find m d ~ 0.05. As in IDutton et al.l (|2007t) . 
we allow md to vary with disk mass in order to fit the 
slope of the RL relation: 



10 103 L. 



0.1 



(8) 



where m^.o and a are free parameters. We tune m^.o and 
a to match the RL relation intercept and slope, respec- 
tively. 

In eq.[6l A' is the effective spin parameter of the system: 



\> = J± X =(J± 



J 2 nnl£| 1/2 



-'200 

GAL 



5/2 



(9) 



2nd 



where jd is the ratio of the disk angular momentum 
Jd to the halo angular momentum J2007 an d A is the 
spin parameter of the halo. Collisionless simulations of 
halo assembly find a log-normal distribution of A that 
is independent of h alo mass. We adop t the distribu- 
tion peak A = 0.042 (jBullock et al.l 120011 ) as our fiducial 
value fo r A 7 . The relaxed ha los in the WMAP5 simula- 
tions of IMaccio et al.l (|2008l ) have a logarithmic scatter 
flog a = 0.228, which we compare to the scatter in A' 
allowed by our RL relation below. 

The parameters / c (c) and /r(A', rrid, c) in eq.|6]account 
for the difference in total energy E between an NFW and 
isothermal halo and the adiabatic contraction of the disk 
in the halo potential, respectively, where c is the halo con- 
centration. We use the exact expression for f c in eq. 23 of 
MMW98, and the approximate expression in their eq. 32 
for Jr. Finally, we compute c as a function of Ld for 
our sample using the best-fitting linear rel ation between 
c and M200 for relaxed, WM AP5 halos in IMaccio et al.l 
(|200l) assuming W = V200/ 2 ijDutton et al.ll2010D . 

Given the different measures of disk size presented in 
this paper, the definition of Rd in cq.|6]merits discussion. 
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The model in cq. [5] assumes that the (cold) baryons in 
disk galaxies are distributed in an exponential stellar disk 
with scale-length Rd- However, it is important to distin- 
guish between Rd, the predicted scale- length of the stel- 
lar mass distribution, and rd, the measured scale- length 
of the projected light distribution. For the variety of 
reasons discussed in §T\ extracting reliable deprojected 
scale- lengths from galaxy photometry is difficult. It is 
therefore not straightforward to relate the measured 
and the predicted Rd- 

For a pure exponential disk with a hxed T^, Rd is di- 
rectly related to any given isophotal radius by the disk 
central surface brightness [1q. If there is no scatter in 
[1q at a given Ld, then the scatter in the distribution of 
isophotal radii is identical to the scatter in Rd- Moreover, 
if the median measured at a given Ld is an unbiased 
estimator of Rd, then the ratio of to the median mea- 
sured isophotal radius relates the latter directly to the 
predicted Rd- 

While the relationship between Rd and isophotal radii 
described above clearly oversimplifies galactic structure, 
the basic scenario meshes with the properties of the 
SFI++. That the scatter in the logi?23.5 — L relation is 
2.5 - 4 times smaller than that previously reported for 
the rd — L relation (see i j3.1|) evidences a disconnect be- 
tween rd and Rd, and suggests that measured -R23.5 are 
better proxies for theoretical Rd than measured scale- 
lengths rd- We embrace this hypothesis in applying eq. [6] 
to the nontemplate RL relation presented in £13.11 

5.2.2. Application to the RL Relation 

We now proceed to tune the parameters of eq. [6] to 
match the RL relation for the nontemplate subsample 
in Fig. [3] We convert i?23.5 to Rd using the best-fitting 
linear relation in Fig. [TJ and convert to solar luminosities 
using Mqj = 4.19. Fig. [T4l shows the RL relation for the 
nontemplate subsample in these units. 

The solid line in Fig. [14] shows eq. [6] with our fiducial 
A' = 0.042 and m d , = 0.054, with a = 0.14. This choice 
of parameters reproduces the slope and zero-point of the 
obs erved relation but is not unique: iDutton et al.l (|2007f ) 
and iGnedin et al.l (]2007l ) demonstrate that a simultane- 
ous fit of the scaling relations, their scatter and their 
residuals is required to break the degeneracies between 
parameters in eq. [6] 

Nonetheless, the implications of the scatter in our RL 
relation for this model can be made clear by consider- 
ing the maximum scatter in each of the parameters of 
cq. [6] that is allowed by the data. We model a con- 
stant scatter e a b s (dashed lines in Fig. [14]) with e mes = 
e RL = 0.04 log(/i- 1 kpc) and e lnt = 0.034 log(/i" 1 kpc) 
(see §4.2|) . The maximum scatter allowed in any one 
parameter of the model in cq. [6] is therefore ej n t . For 
simplicity, we assume that fn is scatter-free in this ex- 
ercise. Eq. IB] shows that the maximum scatters in A' 
and aid allowed by the RL relation are e i og A' = 0-034 
andeiogmd = 0.102, respectively. Since IDutton et al.l 
(|2007[ ) find that the RL relation limits the allowed range 
of these parameters even when simultaneous models to 
all the scaling relations are considered, we consider these 
maximum scatters robust. 

This simple analysis suggests that the range of A' and 
rrid in spiral galaxies is significantly smaller than pre- 
viously estimated: for example, our upper limit on the 
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Fig. 14. — RL relation for the nontemplate subsample expressed 
in scale-length and luminosity units with H = 71.9 km s~ x 
Mpc - . The solid line through the points shows eq. l6l with pa- 
rameters A' = 0.042, rrid = 0.054, a = 0.14. The dashed lines show 
the scatter 2e D (, s = 0.052 log(h — 1 kpc) about that relation allowed 
by the data, and the shaded region shows in comparison the scat- 
ter expected if the distribution of disk spin parameters A' matched 
that of the spin parameters A of collisionless halos. 

sca tter in A' is a factor of ~ 3 smaller than that found 
by IDutton et alJ ([2007), because 6i„t in our RL rela- 
tion scatter is significantly smaller than that of the one 
they model. More importantly, ei og A' is a factor of 6.7 
smaller than ei n y a f ound from cosmological simulations 
(jMaccio et al.ll2008h : to illustrate, the shaded region in 
Fig. [14] shows the expected scatter in the case where 

^log A' — ^log A • 

It is possible that the distribution of specific angular 
momentum jd/md in disk galaxies has less scatter than 
that of their host halos, perhaps due to a redistribution of 
angular mo mentum through bu lge formation. This was 
explored by IDutton et alJ (|2007| ) , but their bulge forma- 
tion scenario reduces the RL scatter by at most 15% for 
realistic bulge-to-disk ratios. It therefore seems unlikely 
that angular momentum redistribution is responsible for 
the narrow range of A' required by our model, particu- 
larly for the Sc-dominated SFI++. 

Another possibility is that disk galaxies form in a 
subset of halos with a distribution of spin parameters 
that is different from that in the halo population as 
a whole. Some studies suggest that this is the case: 
ha los with a quiet me r ger hi story in the simulations 
of iD'Onghia fc Burkertl (|2004h have both lower A and 
ciog a than found for all halos. This possibility has been 
previously discussed in the context of scaling relation 
models (e.g.lde Jong fc Lacevll200Cil: iPizagno et alj|2005t 
IDutton et al.ll2007t IGnedin et alJ 120071 ). However, it is 
unclear whether this effect can explain the extremely 
small range of A implied by our RL relation: once the 
scatter from other parameters such as c and are taken 
into account, the allowed scatter on A' may well be an or- 
der of magnitude smaller than that predicted for A from 
cosmological expectations. 

6. SUMMARY 

Of all the disk galaxy scaling relations, the LV (i.e. 
Tully-Fisher) relation has by far received the most atten- 
tion. It has been used to compute distances to galaxies 
in the nearby universe and study the local velocity field, 
and to put constraints on models of galaxy formation, 
using data at different wavelengths and at both low and 
high rcdshift. While the RL and RV relations can in the- 
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ory serve similar purposes, the large scatter in previously 
published versions of these relations is prohibitive. 

With close to 5000 galaxies, the SFI++ is one of the 
largest available samples to study disk galaxy scaling re- 
lations. The main advantage of the sample is the large 
range of carefully measured parameters available, and 
the well-studied uncertainties associated with them. Ho- 
mogeneous rotation velocities and /—band luminosities 
were published by S07, and the SFI++ LV template re- 
lation studied in detail by M06. In addition to estimates 
of V and L, measures of disk size R are also available for 
a subset of the catalog. It is therefore possible to con- 
struct the full suite of LRV scaling relations from the 
SFI++. 

In this paper, we have presented the LV, RV and 
RL relations for subsets of the template and nontem- 
plate SFI++ samples. Because the SFI++ LV relation 
has been extensively studied, we focus on the RV and 
RL relations. Contrary to previous studies, we adopt 
inclination-corrected isophotal radii rather than disk 
scale- lengths as a measure of the disk size R. While scale- 
lcngth-based studies report an observed sca tter for the 
RL relation of = 0.15-0.2 log(ft. _1 kpc) (jShen et all 
l200l lAvila-Reese et al.ll200l P05, C07), we find e obs = 
0.05 log^/i^kpc). With e obs = 0.11 logO^kpc) and a 
correlation coefficient of r = 0.86, our RV relation is also 
significantly tighter than previously found. We demon- 
strate that the drastically tighter relations and lower 
scatters reported here stem from our use of isophotal 
radii, and argue that they are superior to measured scale- 
lengths because they are less susceptible to measurement 
errors and can be more reliably corrected for internal ex- 
tinction. We examine and interpret the morphological 
dependence of the scaling relations as well as the prop- 
erties of their residuals. 

After carefully accounting for all known sources of 
measurement error and their covariances, we retrieve 
the amount of intrinsic scatter in each of the relations. 
Our error budgets are consistent with constant intrin- 
sic scatter £i nt at high L and high V, with values 
fmt = 0.35 mag, e int = 0.084 log(/i _1 kpc) and e mt = 
0.034 log(/i" 1 kpc) for the LV, RV and RL relations, re- 
spectively. The scatter of the LV and RV relations how- 
ever increases strongly at low V. While this likely reflects 
enhanced intrinsic scatter in V, the effect is difficult to 
quantify due to the competing contributions of peculiar 
velocities, sample incompleteness and measurement er- 
rors. On the other hand, the scatter in the RL relation 
decreases slightly at low L, consistent with the behavior 



of an incompleteness bias. 

We investigate the possible application of the RV re- 
lation as a distance indicator. Comparing with the 
Ccphcid variable distances available for 17 galaxies in 
the template subsample, we find that the RV relation 
returns distances with a median precision of 25%. This 
is slightly inferior to the 15-20% accuracy of LV rela- 
tion distances typically reported. We suggest that RV 
distances are nonetheless accurate enough to provide a 
valid substitute when precise integrated magnitudes are 
not available. We find that the uncertainties on Ho de- 
rived by matching the zero points of the RV relations 
for the template subsample and Ccphcid variable sample 
are comparable to those obtained using the LV relation: 
in both cases, the uncertainties on the Cephcid distance 
zero point dominate. 

To illustrate the potential of the SFI++ scaling rela- 
tions to constrain the formation of Sc galaxies, we apply 
a simple, MMW98-style galaxy formation model to the 
nontemplate RL relation. Assuming that the intrinsic 
scatter in our measured isophotal radii equals that of 
the predicted scale-lengths of the stellar mass distribu- 
tion, we find an upper limit on the scatter of the disk 
spin parameter that is 3 times smaller than previously 
reported, and 6.7 times smaller than the distribution of 
halo spin parameters predicted by cosmological simula- 
tions. As it is unlikely that angular momentum redis- 
tribution due to bulge formation can account for this 
difference, we suggest that the subset of halos in which 
Sc galaxies form has a much smaller spread in spin pa- 
rameters than the broader halo population. A separate 
analysis will be required to quantify the full potential 
of the low-scatter SFI++ RV and RL relations for con- 
straining galaxy formation in light of the improved size 
measurements presented here. 
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APPENDIX 

A. SFI++ SCALING RELATION PARAMETERS AND THEIR UNCERTAINTIES 

In W2.ll wc summarize the derivation of disk luminosities and rotation velocities in the SFI++, and give an overview 
of our approach for computing homogeneous disk sizes for the sample studied here. In this appendix, we present the 
mathematical details of these computations. 

The standard error propagation formula for a function f(x,y) is given by: 

and can be expanded in the expected manner for functions of larger numbers of variables. In eq. IA11 e x and e y are the 
uncertainties on the parameters x and y, respectively, and e xy is the covariance between them. We take covariances 
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into account in the derivation of the SFIH — h error budgets in Appendix B. 

SFI++ Disk Luminosities: Absolute I-Band Magnitudes 

Disk luminosities in the SFI++ are expressed in terms of ab solute J-band magnitude s. As explained in S07 , the 
SFI++ absolute magnitude derivation draws from the work of IGiovanelli et all (|1997bD . lHavnes et al.l (|1999b[) and 
M06; for clarity, we use the notation of S07 where possible. The absolute magnitudes Mi are given by: 

cz 

Mi = m obs - Ai + AM + ki- 5 log -—r{l + z) -25, (A2) 

LlOO/i J 

where m f, s is the observed /-band magnitude extrapolated to 8 disk scale - length s (|Havnes et al.lH999bO . Ai is the 
Galactic extinction ter m from the COBE/DIRBE maps of iSchlegel et al.l (|1998f ). and fcj is a type- and redshift- 
dependent k-correction (|Hanlll992|) : 

ki = (0.1658T - 0.5876)z . (A3) 

We assign distances to the sample galaxi es as described in §2.1.41 The luminosity distance factor (1 + z) in the 5 th 
term on the right-hand side (RHS) of eq. IA21 not included in earlier SFI++ compilations, accounts for cosmological 
dimming. 

The AM term in eq. IA2I corrects for internal extinction, as first derived by IGiovanelli et ail ([1994D : 

AM = 7log(l-e), (A4) 

where e = 1 — b/a is the seeing-corrected ellip ticitv of the disk (IGiovanelli et al.lll997bt lHavnes et al.|[l999bf ) and 7 is 
a magnitude-dependent factor of order unity ([Giovanelli et al.lll995l S07): 

7 = 0.5 for Mi > -19.1 

= 1 - 0.417(il// + 20.3) for - 20.3 < Mi < -19.1 

= 1.0 for - 21.8 < Mi < -20.3 

= 1.35 - 0.35(Mj + 22.8) for - 22.7 < Mi < -21.8 

= 1.3 for Mi < -22.7 (A5) 



We note that there are differences of ~ 0.1 mag between the Mi that we compute using eq. IA2l and the corresponding 
values listed in tables 2 and 4 of S07 (specifically, the values listed in the erratum to that paper). The two main reasons 
for this difference are (1) Mi listed in S07 have the morphological correction of M06 applied. By contrast, we only 
correct Mi for morphology in the LV relation, in the manner described in £13.21 and (2) M06 and S07 a ctua lly use 
a velocity width-dependent value of 7 in computing their Mj, rather than their published relation in eq. IA5I (K. L. 
Masters, priv. comm.). We have verified that these small difference s do not affect the result s presented here. 

The measurement uncertainties e m on Mi are computed as in IGiovanelli et al.l (|1997b| ) . Contributions from cz 
(< 1%) and ki (< 0.01 mag) are small compared to other factors, and are therefore igno red. Assuming that the 
uncertainties e m on to ;, s , on Ai, e 7 on 7, and e e on e arc uncorrelatcd and using eq. IA1[ em is given by: 

4 = 4 + 4 + ^4 + ^^ (A6) 

where 

<J 7 = log(l - e), (A7) 

-7 



iM 



ln(10)(l -e)' 



(A8) 



In the SFI++, e mo b s and e e are computed using the method of lHavnes et al.l (|1999a|) . e 7 = 0.157 ([Giovanell i et al.l 
Il997bf) and e A = 0.2A/ (C. M. Springob, priv. comm.). 

SFI++ Disk Rotation Velocities: Velocity Widths 

Disk rotation velocities in the SFI++ are expressed in terms of velocity widths. As explained in S07, these widths 
are homogeneously derived from either single-dish HI measurements or ORCs, with a preference for the former when 
both are available for a given galaxy. The corrected velocity widths in the SFI++ are given by: 

W =-?-., (A9) 
smi 



where 

Q 



W °»^~ A ° - A f for HI widths, 

l+z 

W bs,ORC 



1+z 



\Y-2 



for ORC widths. 



(A10) 



The HI width measurement technique and corrections are given in ISprineob et all (|2005l) . Briefly, W 2 i, bs 

in eq. IA10I 

is the observed width and A t = 6.5 km s^ 1 is the correction for turbulent motions. The instrumental broadening is 
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A s = 2AuA, where Av is the spectro meter channel separati on in kms -1 and the parameter A depends on the spectrum 

quality and smoothing (see table 2 of iSpringob et al.ll2005D . 

The ORC width derivation is given in iCatinella et al.l (|2005l ) and iCatinella et al.l (120071): as explained in those 
papers, the measured width Wonc,obs in eq. lAlOl is the parametric Polyex model fit (jGiovanelli fc H avncs 2 0021) at the 
radius r opt containing 83% of the total /-band light. Since the HI disks of spirals typically extend to twice this value, 
a correction factor W21/W0RC is applied to homogenize the HI and ORC widths: 

W21 _ f 0.899 + Q.l88r max /r opt for rising ORCs, . . 

Worc ~ \ L075 - 0.013r mQX /r opt for flat ORCs, ^ ll > 

where r max is the measured extent of the ORC. 
The disk inclination i is computed from e using: 

{c0Sl f = ^l_A i ( A12) 

where the intrinsic ax ial ratio qo = 0.13 for galaxies of types Sbc and later, and qo = 0.20 for earlier types 
l|Giovanelli et al.l 11994 S07). Substituting for sini using cq. IA12I and adopting logarithmic units, cq. IA9I can be 
rewritten as: 

\ogW = logQ + 0.51og(l - gg) - 0.51og[l - (1 - e) 2 ]. (A13) 

Assuming that the measurement uncertainties on individual parameters are uncorrelated and applying eq. IA11 the net 
measurement uncertainty ew on log W is therefore: 

4 = el g w = $q£q + %4 + S ew4, (A14) 

where 

*« = M^g- (A15) 



= 1 TTT^T 2v (A16) 



-go 
-(l-e) 



ln(10)(l- [1-e] 2 )' 

The uncertainty cq in Q can be derived from ecis. IA1I and IAIOI assuming uncorrelated errors: 



(A17) 



7ir5^ + 4 t for HI widths ^ 



(1+*) 



u 21 



Wo 1 



for ORC widths, 



where ewobs is the width measurement uncertainty, and the uncertainty on WixIW qrc is assumed to be negligible. 
In the SFI++, e As = 0.25A S , e At = 0.25A t , and e q = 0.15g (jGiovanelli et alj|1997bl) . 

The values of log W and ew computed here are identical to parameters log( Wtf) and e w published in tables 2 and 4 
of S07. 

SFI++ Disk Sizes: Isophotal Radii 

As explained in £12.1.31 we adopt the radius of the disk corresponding to the /ij = 23.5 mag arcsec~ 2 isophote as the 
SFI++ disk size: 

^23.5 = r 23 . 5 [(1-eyP + ARy 1 , (A19) 

w here the terms in bracke ts on the right -hand side correct for internal extinction in the disk following the prescription 
of IGiovanelli et al.l (|1995|) . In eq. IA191 AR is a small perturbative adjustment to account for photometric profile 
variations in galaxies with a given Mj: 

AR = — ln(10) logfl - e)(l - e) a - 2fi ( — - (— >^ , (A20) 

\r23.5 r 2 3.5 J 

(c.f. eq. 7 of IGiovanelli et aT1ll995|) where a and (3 are luminosity-dependent factors shown in their fig. 7. We use the 
best-fitting linear relations: 

a = 0.041 + 1.0M/ , (A21) 

P = 0.031 + 0.863M/ . (A22) 

with Mi from eq. IA2I In cq. IA201 is the measured disk scale-length and (rd/^.s) is the average ratio of disk 
scale-lengths to isophotal radii as a function of Mi for our sample, given by eq. [TJ 
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The parameter r 2 3.5 in eqs. IA191 and [A"20l is the measured isophotal radius r h s .23.5 corrected for galactic extinction 
Ai, a fc-correction kj and cosmological surface brightness dimming per unit frequency interval: 



r 23.5 — r obs, 23. 5 



ln(10)r rf 
2.5 



L4 J -fc / + 2.51og(l + z) 3 ]. 



We use the same valu es of Ai, fej and z as adopted to compute Mj (eq. IA2[) . 
Expressing eg. IA19l in physical, logarithmic units, we obtain: 



logi?23.5 = log(r° 35 ) +log 



cz 



(l0( 



logS, 



(A23) 



(A24) 



where we have adopted the angular size distance and B converts angular size units to radians. We note that the 
cosmological surface brightness dimming term in eq. IA23I partly offsets the cosmological stretching implied by the 
angular size distance in eq. IA24I 

We ass ume that the measureme nt errors on the individual parameters used to derive logi?23.5 are uncorrelated, and 
following I Giovanelli et al.l (|1997bD we neglect the contributions from cz and fcj. We note that because the perturbative 
correction term Ai? is small compared to the measurement errors e r on r {, Si 23.5 for most galaxies in the sample, we 
assume that the uncertainty on this correction simply scales with its amplitude. The measurement errors on logi?23.5 
are therefore: 



2 2 _ s-2 2 

e R — e logR23.5 — °r e r 



fj 2 e 2 



where 



S r = 



>AR e A 



1 



r2 2 

V/3 



°AR e AR "+" °eR e e 



Srd = 



ln(10)r 2 3.5 
A 7 - fcj + 2.51og(l 



$A1 



' 23.5 
7*23.5 



2.5r 23 . 5 



(A25) 

(A26) 

(A27) 
(A28) 

(A29) 
(A30) 

r 23 .5 Vln(10)(l-e) 1 +V (A31) 

Following the approach of iGiovanelli et all (|1997b0 . we adopt uncertainties of ep = 0.15/3 and ear = 0.15A.R. 
Individual measurement uncertainties e r on r t, Sj2 3.5 and e r d on are not available for SFIH — h galaxies. We therefore 
adopt appropriate characteristic values e r = 0.05r o (, Si 23.5 and e r d = 0.15^ (M. P. Haynes, priv. comm.). 



2.5r 23 . 5 
log(l ~ e) 



8ai 



1 



ln(10) 



JeR 



1 23.5 
^23.5 



B. MEASUREMENT UNCERTAINTIES IN THE SFIH 



SCALING RELATIONS 



Below, we derive relations describing the measurement error contributions e mes for the scaling relations constructed 
from R, L and V (see eq. [S]), accounting for correlations between the measurement errors in these quantities. We also 
derive the expected contribution to the intrinsic scatter of each relation from galaxy peculiar velocities. 

The SFI++ LV Relation 

The LV relation in the SFI++ is expressed as a linear correlation between L = Mj and V = log W (jGiovanelli et al.l 
I1997H M06, S07): 

Mi = a LV + b LV log W, (Bl) 

where cllv and b^y are the best-fitting zero point and slope. The net contribution of measurement uncertainties to 
this scatter can be derived from the expressions for Mj (eq. IA2[) and log W (eq. IA13I) : 



<n + 4 + Sfa + {b LV 5 Q Ye 2 Q + (b LV S q ) 2 e z Q + 5, 



J eLV c e 



e ALV 



(B2) 



where S 7 , Sq and 5 q are given in equations IA7| [A15I and [Al6| respectively, and 

$eLV = 5 e M — bLV^eW j (B3) 

where S e M and S e w are given in eqs. IA8I and IA17I (see also IGiovanelli et al.lll997bl M06). The parameter e^LV is the 
uncertainty in the morphological type correction A^y computed in eq. HI which we take as e^LV = 0.15Alv- We note 
that given the number of galaxies in the SFI++, the statistical uncertainties on a^y and &lv are small compared to 
parameter measurement uncertainties, and their contribution to the error budget can be neglected. 
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It is clear that since log W is a distance-independent quantity, peculiar velocities V pec introduce scatter in the LV 
relation only in Mj. Considering the 5th term on the RHS of eg. IA21 a peculiar velocity V pec contributes a scatter 



51og(l + W«0> (B4) 
The SFI++ RV Relation 

We express the SFI++ RV relation as a linear relation between R = log -R23.5 (eq. IA24|) and V = log W (eq. IA13[> : 

log R23.5 = a RV + b RV log W , (B5) 

where a R y and b R y are the best-fitting zero-point and slope. The total measurement error contribution e R y to the 
scatter in the RV relation is therefore: 

e|y = i s = $r4 + S rd4d + + S hl + S ab € ar + (bR.vS Q ) 2 e 2 Q + (b RV S q ) 2 e 2 q + 8 2 eRV e 2 e + e 2 ARV , (B6) 

where 

8 e RV = 8eR — b R vSeW ■ (B7) 

In eq. EH S r , S rd , S AR , 5p, S AR , $q and S q are give n by eqs. IA26I IA271 IA28I IA291 IA301 IA15I and IA16I respectively. 
In eq. IB7[ S eR and 5 e \y are given by eqs. IA31I and IA17I The parameter cahv accounts for the uncertainty in the 
morphological correction of eq. [3l which we take to be 6arv = 0.15A#y. 

As in the LV relation, peculiar velocities introduce scatter in the RV relatio n only along the logi?23.5 axis because 
log W is a distance-independent quantity. The second term on the RHS of eq. IA24I shows that the scatter along this 
axis for non-zero V pec is: 

tRVypec = - l0g(l + Vp ec /cz) , (B8) 

The SFI++ RL relation 
We express the SFI++ RL relation as a linear relation between R = logi?23.5 and L = Mf. 

logi?23.5 = a RL + b RL Mj . (B9) 
The total measurement error contribution to the scatter in the RL relation is therefore: 

4l = d es - %4 + €d4d + tfcl + Sar^ah + &k4 + (b RL 5 7 ) 2 e 2 + S\ RL e\ + 6 2 eRL e 2 e + e 2 ARL , (BIO) 

where 

Sarl — Sar. — b R L , S eR L = S eR ~ b RL 5 e M ■ (Bll) 

In e q. IB101 d r , S r d, 8a r and 5 7 are g iven by eqs. IA261 IA271 lA"29l IA30I and IA7[ respectively, while S eR and S e M 
in eq. IB11I are given by eqs. IA31I and IA8I The parameter €a R l accounts for the uncertainty in the morphological 
correction of eq. [21 which we adopt as carl = 0.15A#£. 

Unlike the RV and LV relations, peculiar vel ociti es in troduce scatter in both variables of the RL relation, since 
log -R23.5 and Mi depend on distance. Eqs. IB8I and IB4I show that a positive peculiar velocity V pec simultaneously 
scatters points towards fainter L and smaller R. The amplitude of the scatter is: 



tRL.ypec = y [51og(l + V pec /cz)f + [log(l + V pec /cz)} 2 = V261og(l + Vpeo/cz) . (B12) 
This scatter vector has a constant slope in the L-R plane: 

~lQg(l + Vpec/cz) _ 

ORL,V P ec - T"j — n TZ -, — r - -U.2 , (B13J 

5l0g(l + Vpec/CZ) 

this value is similar to the best-fitting slope f3 = —0.16 of the measured RL relations (Table [T]). Peculiar velocities 
therefore scatter points almost exactly along the RL relation, and do not contribute significantly to the intrinsic scatter 
in the relation. 



